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We investigate the Mott transition in infinite dimensions in the orbitally degenerate Hubbard 
model. We find that the qualitative features of the Mott transition found in the one band model are 
also present in the orbitally degenerate case. Surprisingly, the quantitative aspects of the transition 
around density one are not very sensitive to orbital degeneracy, justifying the quantitative success 
of the one band model which was previously applied to orbitally degenerate systems. We contrast 
this with quantities that have a sizeable dependence on the orbital degeneracy and comment on the 
role of the intraatomic exchange J. 



ON 

o 



I. INTRODUCTION 

The Mott transition in transition metal oxides has received renewed theoretical and experimental attention. On 
^\ the experimental side, new compounds have been synthesized [1,2] and known compounds such as V2O3 and NiSeS 
t-H have been studied. New insights have been obtained from studying the one band Hubbard model in the limit of large 
Q | lattice coordination [14]. Reviews of the subject have appeared recently [3]. 

1} Several quantitative comparisons between the physics of three dimensional transition metal oxides and the one band 
Hubbard model have been performed by now [4,6] The agreement is surprisingly good, if one considers the fact that 
QO m three dimensional materials such as La,i- x Sr x TiOz and the d 1 electron has a quasi-threefold degenerate level while 
t-H the simplest description of V2O3 involves one electron in a twofold degenerate level [18] . the doping dependence of 
the physical parameters like the effective electron mass in La x Sri- x TiOz can be explained by the one band model 
\Q without adjustable parameters after the values of U and D have been extracted from photoemission data. Similarly, 
r the single band Hubbard model can describe the temperature dependence of both, the optical and the dc-conductivity 
in V2O3 [6]. The photoemission spectra and the optical conductivity of Sri- x CaVOz are also in good agreement with 
the one band Hubbard model [7]. 

In this paper we examine the orbitally degenerate Hubbard model in the limit of infinite dimensions. This publication 
has two goals: a) using the projective self consistent method we show that all the qualitative aspects of the Mott 
transition are independent of the band degeneracy. Furthermore, this method can be used to detect the qualitative 
j3 trends of the dependence of the critical values of the interaction strenght on the band degeneracy, b) To compare 
£^ the predictions of the one band and the two band Hubbard model to justify the success of the one band model in 
explaining experimental data in orbitally degenerate systems. For this purpose we use a recent extensions of the IPT 
(— 1 to particle hole asymetric problems [19]. In the range of densities between zero and one this method is very accurate 
Q and allow us to carry out a detailed comparaison between the one and the two band Hubbard model. We justify a 
O posteriori the success of the one band model for the description of certain quantities related to the low energy features 
of the Hubbard model in the regime of densities < n < 1. We also find that the high energy features of the spectra 
depend on the band degeneracy, as has been known for some time on the basis of less controlled approximations [15]. 

We also explore the dependence of the low energy and the high energy features on the intraorbital exchange 
interaction J which only appears in the orbitally degenerate case. 
Our Hamiltonian for Z& = 2 bands is given by: 

" = ~J= / j c imtr c jm<r H ^~ / j n im<r n ima 
<i j y>,™cr i,racr 

i,m<r i,m<r 

where m = 1,2 (a = ±1) denotes the band (spin) index. A bar over m or a means the complementary value, 
i. e. fa = 2 if m = 1 and m=lifra = 2(<r = —1 if a = 1 and a = 1 if a = —1). 

We introduced three interaction parameters Ui. TJ\ describes the on site interaction between two particles in the 
same band but with opposite spin. TJ2 relates to a pair of particles with opposite band and spin index. TJ3 finally 
concerns the case of equal spin but opposite band index. TJ3 is by an intraatomic exchange energy J smaller than U2 , 
i. e. U3 = U2 — J ■ Rotational symmetry requires U2 = U± — 2J [15]. 
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We study this model on a Bethe lattice with coordination d in the limit of d — > oo. The band density of states in 
this case is a semicircle with a full bandwith 2D = At. 

The lattice model in infinite dimensions is mapped onto an impurity model together with a self consistency condition 

[5,3]. _ ^ ^ ^ ' 

The relevant impurity problem describes an impurity f m(T with band index m and spin a coupled to a bath of 
conduction electrons (cfc ra<7 ). For the frustrated lattice it is given by: 



Himp — €f y ] fmafma + ^ ] e k c ^ m(T c km(T 



y "j { c kmaf m < T fma c kmtr) + ^ ^ ] n ml n mcr 
k,mcr ma 



U l V n if) + — V n (/ } (2) 
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where the hybridization function A(w) := L0 _ e *_ irl has to fulfill the self consistency condition: 

A(w)=t 2 G(w) (3) 

The fact that the general structure of the model is very similar to the one band case, can be understood on the 
basis of simple physical arguments along the lines which have already been presented in the one band case [3] and the 
projective self consistent method. We discuss this in section II. To obtain quantitative results we use a generalization 
of iterative perturbation theory (IPT) to particle hole asymmetric problems which we introduced earlier [19]. This 
approximation is more accurate than 10 % in the one band Hubbard model, and captures all the qualitative features. 
The scheme is based on an interpolating self energy which becomes exact in various limits, namely for low and high 
frequencies as well as in the atomic and the small U limit. The comparaison between the one band and the two band 
results based on the IPT is presented in section III. 

There have been several recent investigations of the degenerate Hubbard model using a variety of techniques, the 
Gutzwiller approximation [22], the slave boson technique [21] and the quantum Monte Carlo method [23]. On the 
qualitative level our results are in agreement with these studies. 



II. THE MOTT TRANSITION IN A DEGENERATE HUBBARD MODEL: PROJECTIVE SELF 

CONSISTENT ANALYSIS 

The degenerate Hubbard model is insulating at integer particle number per lattice site provided that U is sufficiently 
large. Figure 1 displays the total particle density per lattice site as a function of the chemical potential for {7 = 4 
obtained with using the IPT method. There is a Mott transition at n tot = 1, 2, and 3. Notice that the gap around 
the ntot = 1 transition is larger than the gap around the n tot = 2 transition. In addition, we found the critical value 
U c for the instability of the metallic solution is larger for n tot = 2 than for n tot = 1. The transitions at n tot = 1 and 
ntot = 3 are equivalent because of particle hole symmetry. 
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FIG. 1. Particle density of the two band Hubbard model as a function of fi for U = 4. 
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In the one band Hubbard model, in infinite dimensions, the nature of the Mott transition remained a controversial 
issue for some time and was finally resolved by the application of the projective self consistent method. [9] [16] We 
know now that at zero temperature there is a region in the chemical potential and interaction plane where a metallic 
and an insulating solution coexist. The zero temperature metal to insulator transition is second order and is driven 
by the divergence of the effective mass. 

We have extended the projective self consistent method to the degenerate case and we show that the nature of the 
Mott transition is not qualitatively changed by orbital degeneracy. 

Simple approximations within this scheme explain the differences in the critical U c and in the size of the band gaps 
that we observed in IPT and which have also been found in QMC calculations [23]. 

In the projective self consistent scheme one separates the set {k} of "bath" electrons into sets H C {k} and L C {k} 
with corresponding couplings {ej? , V k H } and {e k , V k L } describing the high and low energy features. Then one projects 
out the high energy degrees of freedom to obtain an effective self consistent problem for the spectral density, involving 
the low energy variables only. 

The weight of the quasiparticle resonance is w = YlkeL (^k) • ^or notational clarity we also affix an L to e k and 
c ka for k G L. The quantity w is the small parameter in a singular perturbation theory. 

To construct the exact conditions for the Mott transition we separate the impurity Hamiltonian into three parts as 

~Ham = ~H a +~Hb + 7i m (4) 
The high energy ("atomic" -like) part H a is given by 

' Ha = \ E n m*n m , a , + £ V k H (fa c ka + h - c -) 

(m<T)^(m> <t>) a,keH 

+ £ (<* ~ ri( C ka) + Cka (5) 
a,keH 

H a is the Hamiltonian of the Mott insulator with proper occupancy. It has d degenerate ground states carrying a 
representation of SU(N) (d = 21}, for n tot = 1 and d = Z5 — 1) for n tot = 2; N = 21^). We label the ground states 
as \l> a and define Hubbard operators = \l > a a < l'\ acting on this low energy manifold. 

The low energy spectrum of the impurity Hamiltonian consists of this degenerate manifold combined with excitations 
of the low energy bath electrons described by the Hamiltonian 

Tib = £ e k L (c L ka y c L ka (6) 

a,k£L 

These low energy states are mixed with the high energy excitations of 7i a by 

Km = £ V k L (f+ct + he) = V LH + V HL (7) 

a,k£L 

since f\ or f a acting on \l > a create particle or hole states with energies increased by at least the particle or hole 
gap. It is convenient to define normalized operators of the local low energy bath electrons at the impurity site, 

keL 

The total Hamiltonian is then of the form Ham = ( -i,hl u I • Applying a canonical transformation, one obtains 

V " H a J 

an effective Hamiltonian which is block diagonal and does not mix low and high energy states: ~H e ff '■= g s H.amg~ s = 
"H^ \ 

eff h )■ Here, the effective Hamiltonian Ti.^ s is defined on the low energy Hilbert space {\l > a } ® {k £ L}. 



n& 

j 

F£ L = V^X lv J?« Cfl where 



Analogously, the operators f a are transformed into F a := e f a e = I j^Jil jpHH )• ^ n particular, we obtain 

\ a a 

-m w -V- (9) 
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and summation over repeated indices is assumed 

The indices /i, /i', and a in (9) denote the 2Z& spin orbital combinations (ma). The same is true for all other indices 
as far as n tot = 1 is concerned. If n tot = 2 and Z& = 2, however, the indices Z, Z' , p, and p' should be thought of as 
double indices, e.g. Z = (Zi, Z2), describing one of six possible pairs (ma, m'a'), which denote the quantum numbers of 
the ground state. 

If ntot = 1 (h = 1 or Zj = 2), the J^f are of the form J^f = ai Su'S^a + bi SvaSi^, while for n tot = 2 (Z5 = 2) one 
has J/;, a = a 2 VA'^« - b 2 eP al ^eP^. 



The low energy part of the self consistency condition becomes t 2 G LL (u>) = X^fcgi v-J L ' wnere G LL (t — t') 



TF LL {t) {F LL {t r )) + \ 

In addition, one can write the effective low energy Hamiltonian as 

n^ = n b + wD 2 [J s O s +J p O p }+ const (10) 

where O s and O p are the SU(N) (N = 2Z&) generalizations of the spin spin interaction and the potential scattering 
between the impurity states and the low energy part of the bath of conduction electrons, respectively. The parameters 
J s , J p are linear combinations of the exchange matrix elements in eq. 9 and will be discussed below. For n tot = 1 and 
SU(N = 2l b ), O s has the form 

O s = {x illl ,c ll c+ - ^X^c^c+^J (11) 
which for Z;, = 1 (i. e. TV" = 2) can be written as O s = 4 Sc^a^^ic^i . For lb = 2 and n tot = 2, O s is given by 

O s = (^X pp c,c+ - e^e^'^X pp , c„c+ ) (12) 
The potential scattering is described by 

O p = X^^c+Cf, (13) 

In equations (11) to (13) as well as below summation over repeated indices is assumed. 

The condition for the destruction of the metallic solution follows from the self consistency condition and is given 
by t 2 ({(F^ l ) + , F% L }^ = Y*k£L ( V k L ) 2 - Evaluating this expression we obtain: 

^ = 2Zft ^PP lC ^ c t''} {^p'' (yJp'i' ) ~ 'Kp< (jip ) } 

+ h £ «• |2 ,14) 

W ficx 

It assumes a very simple form in specific cases: For n to t — 1 and ^ = 1 

± = (l-( 0s ))J 2 +J p 2 (15) 
with J s = bi and J p = ai + for n to t = 1 and h = 2 



J_ _ 15 
t 2 ~ ^16 



2-(H-(°s)) J s+ J P 2 ( 16 ) 



with J s = &i and J p = ai + &i/4, and for n tot = 2 and lb = 2 

^=(1-(°°)) J s+ J2 P ( 17 ) 

with J s = &2 and J p = a,2 — ^f. Here (...) indicates the expectation value with respect to the effective low energy 
Hamiltonian Ti.^ s . 
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An expression for the parameters a 8 - and 6 8 - (and therefore J s and J p ) can be obtained by an expansion in the 
hybridization. It turns out that the main differences between the transitions at n to t = 1 and n to t = 2 appear to 
leading order in this expansion. To illustrate this point let us evaluate J^f for the single band model at n to t = 1 
and for the two band model at n to t = 1 as well as n to t = 2 (and J = 0). We obtain for n to t = 1 (lb = 1 or Z& = 2) 
ai = — Kv and 6i = — + TT X ■ For n* n * = 2 and h = 2, in contrast, one has an = — Kr and 6o = — + tft — • These 
results can be used to estimate U c i, the critical value of the interaction where the insulator ceases to exist, as well as 
the regime of the stability of the insulating solution in the /i — U plane. As discussed in reference [16] these quantities 
are obtained by noticing that for small perturbations (Xppic^c^,) ps (X pp i) {cpC^, ). The resulting conditions can be 
used to determine for given value U the size g\ of the regime where the insulating solution is stable ("insulating gap"). 
Around particle hole symmetry, i. e. for n to t = lb this quantity is given by g\ = 2A/i c , where an analytic expression 
can be obtained for A/i Cl : 

a "" = (f) S + ¥( 1 -\/ 1 + 2 « fc + 1 »(f) S ) « 18 » 

We see that the insulating gap in the two band model at n to t = 2 is smaller than that of the one band model at 
n to t = 1. In addition, it can be shown by evaluating the equations numerically that the gap for lb = 2 and n to t = 2 is 
also smaller than that for lb = 2 and n to t = 1. A plot of the critical /i Cl as a function of U is shown in figure 2. When 
comparing the insulating gaps for n to t = 1, we find that the values are very similar for lb = 1 and lb = 2. 

For the critical interaction at which the insulating gap closes we obtain U Cl (lb = 2,n tot = 1) = 1.63 < 
U Cl (lb = 1, n tot = l) = V3<U Cl (lb = 2, n tot = 2) = VE. 

Notice that U Cl (lb = 2,n tot = 1) = 1.63 is only slightly below U Cl (lb = l,«tot = 1) ~ 1.73. 




In the large d Hubbard model, the single particle gap or more precisely the jump of the chemical potential when 
going from hole doping to electron doping is smaller than the optical gap in the insulator [16]. We will denote this 
gap by 32 and it can be extracted from the n vs [i. 

A numerical determination of 32 requires the evaluation of the expectation values (O s ). Experience with the one 
band case has shown that these expectation values are in the intermediate coupling range. A detailed numerical 
evaluation of these is left for a future publication. It is instructive however to understand the trends that appear as 
we change band dgeneracy, by taking one half of their values in the strong coupling (sc) limit, an assumption which is 
very accurate in the one band case [9] [16] . The strong coupling limit can be calculated analtyically < O s > sc = — § 
for ntot = 1 and l b = 1, < O s > sc = —5/4 for n to t = 1 and lb = 1, as well as < O s > sc = —5 for n to t = 2 and lb = 2. 
To estimate 32 as well as the critical value U c . 2 at which the metallic solution breaks down, we assume therefore 
< O s >= I < O s > sc - This yields U C2 = 2.4 for lb = 1 and n to t = 1> U C2 = 2.5 for lb = 2 and n to t = 1> and U C2 = 3.9 
for lb = 2 and n to t = 2. Notice, that the values for n to t = 1 are again very similar. To decide which value is smaller in 
the exact solution, the computation of the first correction in jj along the same lines as in the one band model [17] are 
necessary. The critical interaction for lb = 2 and n to t = 2, in contrast, is clearly larger than for n to t = 1- The region 
in the [i — U plane where the metallic solution is stable is indicated in figure 3, where we plotted A/i C2 as a function 
of U. For n to t = 1> this quantitiy is defined by A/j C2 := fi C2 — \U while for n to t = 2, A/j C2 := fi C2 — | U. 32 depends 
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only weakly on orbital degeneracy, for {7 = 4 and n tot = 1 we obtain, for example, g 2 = 2.6 for Z& = 1 and g 2 = 2.4 
for lb = 2. On the other hand, in the two band model this gap is found to be smaller at n tot = 2 than at n tot = 1 
(g 2 = 0.7 for l b = 2, n tot = 2, and U = 4). 
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FIG. 3. Stability of the metallic solution: A/i C2 as a function of U for a different number of bands h and fillings ntot = 1 or 
2. For ntot = 1, A/t C2 := /t C2 — \JJ while for n tot = 2, A/t C2 := /t C2 — |Z7. 

All the qualitative results obtained for {7 C; and gf, are in agreement with IPT calculations discussed in the next 
section and with the recent exact QMC calculations of ref. [23]. 



III. IPT: COMPARISON BETWEEN THE ONE BAND AND THE TWO BAND MODEL 

The Mott transition as a function of filling (from the metallic side) is driven by the collapse of an energy scale, the 
renormalized Fermi energy and the consequent divergence of the effective mass. The control parameter is either the 
interaction U or the chemical potential [i. The goal of this section is to compare the results of the two band against 
the one band model illustrating that many quantities depend weakly on degeneracy, which justifies the success of the 
one band model in comparing it with experiments in transition metal oxides. We also notice that the high energy 
features of the spectral function have a sizeable dependence on band degeneracy in the range of densities between 
zero and one. 

The comparaison is done using the recent extension of the IPT scheme to particle hole asymmetric problems. This 
extension is very accurate for denisties between and 1. 

The evolution of the spectral function of the two band model with increasing filling is displayed in figure 4 for 
ntot < 1- We observe a quasiparticle peak whose width determines the quasiparticle residue Z and scales to zero as 
the Mott transition is approached. 
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FIG. 4. Evolution of the spectral function of the frustrated two band Hubbard model with increasing filling for ntot < 1 and 
[7 = 4(J=0). 

Figure 5 shows the doping dependence of Z. We added the one band result for comparison. The data for the two 
band case is very close (and slightly smaller) than those for the single band model. The quasiparticle residue vanishes 
as the Mott transition is approached indicating a divergence of the effective mass. 
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FIG. 5. Quasiparticle residue Z = as a function of the total particle number ntot for U = 4 and J = 0. 

In figure 6 we show the {/-dependence of the quasiparticle residue at finite doping. It decreases from its noninter- 
acting value (i. e. Z(U = 0) = 1) but remains finite for large {7. 
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FIG. 6. Dependence of the quasiparticle residue Z = on the interaction U for ntot = 0.8 and J = 0. 
At n tot = 1, in contrast, Z vanishes above U c (see figure 7). 
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FIG. 7. Comparison of the Z vs U plots for the two band and the single band model at ntot = 1 (J = 0). 



In comparison with the single band model, the slope of the n vs [i plots is steeper for the case of two bands (figure 

1 dn 
n 2 dfi 



8). Thus the compressibility re = §77 is increased as a result of orbital degeneracy. 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



U=4 




-1.0 -0.5 0.0 0.5 



1.0 



FIG. 8. Comparison of the tivs(i plots for the two band and the single band model (U = 4, J = 0). 

It is clear that the low energy parameters are not a strong function of the degeneracy as is seen, e. g., from the 
similarity in the Z vs U curves shown in figures 6 and 7. This is different from the high energy features. Notice 
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the differences in the weights of the Hubbard bands (figure 9) which should be readily observed in photoemission 
experiments. At n tot = 1 the weights of both Hubbard bands are equal (i. e. |) in the single band model, while in 
the two band problem the weights are ^ and | respectively. 
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FIG. 9. Comparison of the spectral functions in the one and the two band model for U = 4: a) Mott insulator (ntot = 1)- 
b) metallic state for small doping (ntot = 0.94). The chemical potential for the Mott insulator in figure a) is chosen such that 
the difference fiiband — fJ-2band is the same as for the metal plotted in figure b. 

A more precise probe is given by the optical spectroscopy as discussed in the following section. 

IV. OPTICAL CONDUCTIVITY 

Figure 10 contains the optical conductivity <t(u>) corresponding to the plots of figure 9. The data were calculated 
from: 

<r(w)=l b <T / dep a (e) I —A e (w') 

U J-oo J-oo ^ 

A e (w' + w)(n f (w') - n f (w' + w)) (19) 

where <t = and A ek (w) = — 2ImGt(w) = -2Im ^ +)[ _ (; ^ s <(^) ' parameter li denotes again the number of 
bands (Z& = 1 for one band and Z& = 2 for the two band model). 

For the insulator (figure 10 a) we find a single maximum around u> = {7. The optical conductivity of the metal 
at small doping (figure 10 b), in contrast, is more complicated. In both, the one and the two band case, there are 
three main maxima Pi, P2, and P3. Pi stems from excitations from the lower Hubbard band to the resonance peak, 
P2 represents excitations from the resonance peak to the upper Hubbard band, and P3 corresponds to excitations 
from the lower to the upper Hubbard band. This peak is split into two maxima due to the satellite peaks in the 
lower Hubbard band. The smaller feature between the Drude peak and Pi for Z& = 2 is due to the satellite peaks 
which appear in the spectral function next to the resonance. In the two band model these satellite peaks are quite 
pronounced. Notice that Pi has almost the same weight in the one and the two band model because the factor 2 for the 
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orbital degeneracy compensates for the fact that the weight of the lower Hubbard band is reduced by approximately 
|. At sufficiently high temperatures all the fine structure disappears and our curves for the one band model are very 
close to the previous Montecarlo studies by Jarrell et. al. [20]. 
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FIG. 10. Comparison of the optical conductivity of the one and the two band model for U = 4 (and J = 0): a) Mott insulator 
(ntot = 1) and b) metal at small doping (ntot = 0.94). 

In figure 11 we compare the optical conductivity of the one and two band Hubbard model with one electron per 
site at U = 2. The intensity of the mid infrared peak differs since it originates from both, excitations between the 
resonance and the lower as well as the upper Hubbard band. The qualitative features however are not changed by 
orbital degeneracy. 
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FIG. 11. Comparison of the optical conductivity for the one and the two band model at U = 2 and ntot = 1. 
The optical conductivity of the model, at T = can be represented as: 

uj* 2 

<t(w) = -p-c5(w) + c» e <,H (20) 

47T 

where the coefficient in front of the delta function is the Drude weight and u>* p is the renormalized plasma frequency. 
The Drude weight and the integrated optical spectral weight evolve as a function of doping. In infinite dimensions the 
Drude weig ht is simply related to Z [24,25] via: ^ = 4?rt n %J h Zp°(0) The integrated optical weight is proportional 
to the kinetic energy < T >. We evaluated the contribution from <j reg (u>) for the cases shown in figures 10 b) and 11. 
We find that the total spectral weight at finite frequencies w reg := J Q a reg (u)du is larger for the two band model 
than in the single band case: "W'*^) ~ 1.45 f or jj = 4 an d n tot = 0.94 (figure 10 b) and "W'*^) _ Lg4 im jj = 2 
and ntot = 1 (figure 11. 

The evolution of the optical conductivity of the doped metal ({7 = 4) with increasing filling is displayed in figure 
12. The data shown correspond to the spectral functions of figure 4. For the dilute system <t{u>) does not have much 
structure. It increases for w — > and has a small maximum at high frequencies corresponding to excitations into the 
upper Hubbard band. As the filling is increased, the resonance peak develops in the density of states. This results 
in the formation of the associated peaks Pi and P2, indicating excitations between the resonance and the Hubbard 
bands. At the same time the weight in the upper Hubbard band increases. 
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FIG. 12. Evolution of the optical conductivity with increasing filling (U = 4, J = 0, h = 2). 

Katsufuji et al. [12] obtained similar plots for a whole series of substances Ri- x Ca x TiOz with rare earth ions 
R = La, Pr, Nd, Sm, and Y. They observed that an isosbectic point u> c where the optical conductivity vs doping 
curves intersect, and computed tn effective number of carriers: 

N *ff ■■= -Z^T / <") d " ( 21 ) 
we Jo 

Fitting a linear behavior to the experimental data for small doping, they showed that the corresponding slope C 
increases drastically as U c is approached: C ~ 1/ (U — U c ). 

Motivated by these experiments we computed the same quantitiy in the degenerate Hubbard model, for U = 
2.5, 3, and 4 (figure 13). In all three cases the curves are not really linear. The mean slope is enhanced near U c = 2.9. 
However, comparing the data for {7 = 4 and U = 3, we see that the slope does change very much as U c is approached. 
These findings are in agreement with Monte Carlo data for the two dimensional Hubbard model [13]. 
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FIG. 13. Neff/a 2 for different U. The calculations were done at T = 0.02 with a finite T = 0.2. 

V. EFFECTS OF A FINITE INTRAATOMIC EXCHANGE J 

An important parameter with no analogue in the one band model is the intraorbital exchange J, which originates in 
Hund's rule. This parameter does not affect the low energy features significantly but affects the high energy properties 
as shown in figure 14. Once J is finite, one is dealing with three different interaction constants (Ui = U , U2 = U — 2J 
and U3 = U — 3 J) causing a splitting of the upper Hubbard band. 
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FIG. 14. Evolution of the spectral function for U = 4 and J = 0.5 with increasing particle density ntot- 
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The interatomic exchange has a profound influence on the phase diagram. A detailed calculation of the competition 
between the different possible phases depends very sensitively on the nature of the lattice and will not be undertaken 
here. Instead we will expand on an earlier observation [19] that the one of the main roles of the orbital degeneracy 
in the context of the dynamical mean field theory is to stabilize over a wide temperature range the paramagnetic 
solution. In the two band Hubbard model the energy difference between the the ferromagnetic state and the Neel 

t' 2 T 

state is given by the energy scale jjt, which can be rather small since J tends to be much smaller than U. Since 
the difference between the various ordered configurations is smaller in the orbitally degenerate situation the system 
remains in the paramagnetic phase over a wider range of temeperatures. 

We present below a simple strong coupling argument as to why this this effect is more pronounced in realistic 
models which incorporate the hopping of the transition metal ion via the ligand hole. A more quantitative study 
of the ordered phases of the Y\- x La x TiO?, system has been carried out recently by Mizoaka and Fujimori using an 
expansion around the unrestricted Hartree Fock approximation. [26] 

We start with a model of LaTiOz containing three Ti 3d orbitals degenerate orbitals (xy, xz, yz). hopping via a 
filled oxygen shell at an energy (e d — A) . The complete Hamiltonian is the sum of three terms: 

H site = e d d La d im<r + (e d " A) Pt+6,l,a Pi+SM ( 22 ) 

Here, i labels the position of the Ti sites. The index 6 is half a unit vector 

(S G {^e x , ^e y , ^e z }) so that i + S marks the position of the oxygens. Moreover, m and I are band indexes: 
m G {xy, xz, yz} and I G {x, y, z}. 

There are only two pairs of orbitals which overlap: the Ti Zd xy with the O 2p y and the Ti Zd xz with the O 2p z (xz 
plane). This results in the kinetic energy: 

H ki „ = t pd ( d txy,a Pi+^e^y,* + h.C. + d+ xy(T p i+ l ey , X)(T + h.C. + 

«>;7G{-1,1} 

d t,xz,a Pi+%e x ,z,a + h.C. + df xzo . Vi^ 2 e z ,x,a + h.C. + 

d t,yz,a Pi+^e y ,z,a + h.C. + df yz (T p i+ X e ^ y ^ + h.C.) (23) 

We model the on site interactions by a degenerate Hubbard model: 

Hint = U± ^ niraf«,iraX + U2 ^ WiratWira'X + U3 ^ Wiracr«.ira'cr (24) 

Here, Ui is the on site interaction at a Ti site between electrons in the same band and opposite spins. Analogously, U2 
(U3) couples electrons in different bands and with opposite (parallel) spins. The parameter U3 is by an intraatomic 
exchange energy J smaller than U2 (U3 = U2 — J, J > 0). Moreover, rotational symmetry requires Ui — U2 = 2J. We 
neglected the on site interaction on oxygen sites as well as next neighbor interactions. The sums over the spin and 
band indices are to be understood such that each pair of states (ma, m'a') or (la, I' a') is included only once. 

H = H s it e + H]~in + Hi n t (25) 

We derived the effective Hamiltonian describing the spin excitations of the insulator by using perturbation theory up 
to fourth order in the O-Ti hopping t pd . The procedure is tedious but straigthforward and is most easily summarized 
by giving the energies of the different magnetic configurations of a pair of neighboring Ti atoms as listed in figure 15. 
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FIG. 15. Independent configurations yielding different energies. 
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FIG. 16. Energy per bond for different configurations. (The configurations are shown in figure 15). 

We used a realistic set of parameters (Ui = ieV, A = 6eV and t pc i = O.beV). The qualitative results, however, 
are not very sensitive on the actual choice of the numerical values. Notice that the energy difference between the 
lowest energy configurations extremely small (less than 5 mev). We checked that in the limit of extremely large A 
the energy gets minimal for ferromagnetic bonds for all values of J. In this case the ferromagnetism of a two band 
Hubbard model is recovered. 

Clearly over a large range of the intraatomic exchange J an almost non magnetic configuration, i. e. a configu- 
ration with a very small magnetic exchange constant, is lowest in energy. For extremely large J the ferromagnetic 
configurations are favored while the antiferromagnetic contributions related to Ji have always higher energies. For 
LaTiOz, we conclude from our simple model that antiferromagnetism is strongly suppressed and the paramagnetic 
phase is stabilized justifying our use of the paramagnetic solution in our earlier descrciption of the experimental data 
in transition metal oxides using the one band Hubbard model. 

It is clear that magnetic ordering takes place only below a quite low transition temperature T c . Above T c the 
system is more or less frustrated. In this temperature regime the dynamical mean field theory has to be applied to a 
fully frustrated lattice to yield a reasonable description. 
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